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I. INTRODUCTION 


The prospect |IJ.|2| of observable baryon number violation in high energy elec- 
troweak collisions has provoked much excitement, but despite considerable effort, it 
has not been possible to obtain conclusive evidence that unsuppressed baryon num¬ 
ber violation can indeed occur at energies in the multi-TeV range. The purpose of 
this paper is to summarize some recent developments which offer hope for a reliable 
calculation of the magnitude of baryon number violating cross sections. 

Limitations of space do not permit a comprehensive treatment of this vast subject. 
These lectures will focus on the semiclassical methods introduced in Refs. [3] and [|] 
and further work by the present authors. The main idea in Refs. 0 and 0 consists 
in studying processes with energy and initial particle number of the form E = e/ g 2 
and Ni = z//g 2 , with g being the electroweak coupling constant, in the limit where 
g —> 0 while e and v are held fixed. Under such conditions one can justify the use of 
semiclassical methods to extract the exponential behavior of the semi-exclusive cross 
section, which can be done by calculating hybrid Euclidean-Minkowski solutions which 
violate baryon number via tunneling from one (non-vacuum) real-time configuration 
to another. Finding such such semiclassical solutions is however highly non-trivial and 
can only be done by computational methods. In Ref. 0 we have applied numerical 
techniques to the study of processes in which the gauge and Higgs fields change 
their topology through purely classical evolution, thus inducing a violation of baryon 
number. The investigation of the semiclassical solutions that account for tunneling 
under the barrier is in progress. 

These two approaches, one involving classically forbidden tunneling-like processes 
and the other classically allowed topological transitions, probe complementary aspects 
of the problem and should produce compatible results (thereby providing consistency 
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checks). Both approaches require solutions to certain nonlinear partial differential 
equations for which no (nontrivial) exact solutions are known; however, these equa¬ 
tions are well suited to computational study and one can still make considerable 
progress. Along with P. Tinyakov, we have recently launched a numerical investi¬ 
gation of the standard model baryon number violating rates based on Ref. [[J, and 
in the next section we shall review the relevant formalism and present some of our 
initial results. The rest of this paper is devoted to the numerical work of Ref. ||. 
We presently have much more to say about this approach because its correspond¬ 
ing computational study is at a mature stage of development. Finally, in an effort 
to write a self contained work, the remainder of this introduction is devoted to a 
brief exposition of nonperturbative baryon number violation in the standard model 
of electroweak interactions. 

For our purposes, when we talk of the “standard model” we mean the standard 
model in which the Weinberg angle has been set to zero, i.e. we shall be considering 
577(2) gauge theory spontaneously broken via a single Higgs doublet. This simplified 
model has all the relevant physics. Most important, the gauge structure dictates non¬ 
trivial topology for the bosonic vacuum sector. Working in the temporal gauge with 
periodic boundary conditions at spatial infinity, each vacuum may be characterized 
by an integer called the winding number which measures the number of times the 
gauge manifold is wound around 3-space ||. As this number is a topological invari¬ 
ant, vacua of different winding numbers cannot be continuously deformed into one 
another. 

Because of the axial vector anomaly, baryon number violation occurs when the 
gauge and Higgs fields change their topology 0. Adjacent topological sectors are 
separated by an extremely high barrier, the top of which is a static saddle-point 
solution to the equations of motion. This configuration is called the sphaleron 0, 
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and it has an energy of about 10 TeV and possesses a single unstable direction in field 
space. At low energy the baryon number violating rates are exceedingly small, as 
the gauge and Higgs fields must first tunnel through the sphaleron, which is indeed 
extremely unlikely. 

Recently, the prospect of rapid baryon number violation at high temperatures 
and high energies has emerged. The basic idea is that if the gauge and Higgs fields 
have enough energy, they can change their topology by sailing over the sphaleron 
barrier rather than being forced to tunnel through it. At high temperatures this 
is precisely what happens, and it is generally agreed that baryon number violation 
becomes unsuppressed in the early universe ||. 

The situation in high energy collisions is far less clear. The limiting process in 
baryon number violation is the production of a sphaleron-like configuration from an 
incident beam of high energy particles. But since the sphaleron is a large extended 
object, there is a scale mismatch with the initial high energy two-particle state, and 
hence one naively expects the baryon number violating rate to be small. However, 
Ringwald 0 and Espinosa || have suggested that the sum over many-particlc final 
states gives rise to factors that grow with energy sufficiently rapidly to offset any 
exponential suppression. If true, this offers the exciting prospect of one day observing 
baryon number violation in high energy collisions. 

II. THE CLASSICALLY FORBIDDEN DOMAIN 


The approach of Ringwald and Espinosa [l]||], however, neglects some important 
corrections which still elude calculation despite considerable effort. Apart from lattice 
simulations, semi-classical techniques are our only handle on nonperturbative effects. 
The basic problem with calculating the anomalous baryon number violating cross 
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sections is that exclusive two-particle initial states are not very amenable to these 
methods, and there are potentially large initial state corrections whose effects remain 
undetermined. 


Rather than calculating the two-particle cross section directly, Rubakov, Son and 
Tinyakov jdj] investigate a related quantity for which semiclassical methods are still 
applicable. Their method involves saturating the path integral representation of this 
quantity with a complexified Euclidean-Minkowski saddle-point. This solution in¬ 
cludes the effects of tunneling under the sphaleron barrier, and is a generalization of 
the periodic instanton of Ref. @. We now review in more detail the work of Ref. Q, 
along with selected portions of Refs. H - |12| upon which this work is based. 


A. The Inclusive Cross Section 


As previously mentioned, the calculations of Refs. |l| and @ for the two-particle 
baryon number violating cross section, cr 2 (P), become unreliable at high energy. This 
is because of a failure of semiclassical methods in calculating exclusive quantities like 
two-particle scattering amplitudes. So rather than calculating ( 72 (E) directly, Ref. 
examines a related inclusive quantity: 

a(E,N) = J2 \<f\SP E P N \i>\ 2 , (2.1) 

/,* 

where the sum is over all initial and final states, 5 is the 5-matrix, and P E and Pjy are 
projection operators onto subspaces of energy E and particle number N respectively. 

Unlike the exclusive two-particle cross section, <j(E, N ) is directly calculable by 
semiclassical methods as long as the incident particle number N remains large. If the 
energy and particle number are parameterized by 

E = 4 (2.2a) 

r 

N = ^ , (2.2b) 

9 
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and the fields are rescaled by appropriate powers of g 2 , then in the limit g —■► 0 with e 
and v held fixed, the path integral for cr(E, N ) can be saturated by a ^-independent 
classical saddle-point solution to the equations of motion. As shown in the next 
section, the cross section takes the form 


cr(E, N ) = exp 


F(e, 0 + 0 ( 9 °) 


(2.3) 


where the function F(e, v) is determined by the classical solution. 

The utility of a(E } N ) is that it may be used to bound the two-particle cross 
section and allow one to extract the exponential behavior of a 2 (E). By construc¬ 
tion, cr(E,N ) provides an upper bound to a 2 (E). This is because one of the initial 
iV-particle states of ( |2.1|) possesses N — 2 free propagating particles and two collid¬ 
ing particles [[l]]]. A lower bound may be obtained under some reasonable physical 
assumptions |L3|. Let \^>n > be the initial state that saturates the sum in ( |2.1|) . If 
the process 2 —> any proceeds through some preferred intermediate state, such as a 
sphaleron-like configuration in the case of baryon number violation, then the sub¬ 
stitution of this state by \ipN > will underestimate the result, giving the inequality 
|< 'ipN 12 >| 2 a(E,N) < cr 2 (E). Estimating |< ip N |2 >| 2 ~ exp (—const N), together 
with the previous upper bound, gives the inequalities 


exp (—constN)a(E,N) < a 2 (E) <a(E,N) , (2.4) 

from which it follows that 


lim g 2 lncr 2 (E) = F(e, u) + 0(u) . (2.5) 

The consistency of the first inequality requires that F(e, v) has a smooth v —> 0 
limit, in which case F(e, 0) determines the exponential behavior of a 2 (E). However, 
the second inequality of (|2.4|) holds regardless of continuity, and hence if cr(E } N) is 
exponentially suppressed (for any value of N), then so is <j 2 (E). 
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B. Development of the Formalism 


We now review the formalism developed in Ref. Q used to calculate the inclusive- 
state cross section a(E,N). For purposes of illustration, we consider a system 
with a single real scalar field, whose field operator 0(x) has eigenstates defined by 
0(x) |0 >= 0(x) |0 >. The approach that follows is based on a coherent state for¬ 
malism, where coherent states |a > are defined by a k |a >= a k |a >, with a k being 
the annihilation operator of the k-th mode. In field space the coherent states take 
the form 


< 0 |a > = const ■ exp 


a k 


i ' ^ / i \ 


where 0(k) is the spatial Fourier transform of 0(x), given by 

^ (x) = /(^ e ‘ k ^ (k) ' (27) 

In the coherent state formalism, the 5-matrix is represented by its kernel 
S(b*,a) =< b\S\a >, and inserting a complete set of field-states we can write 


S(b*,a ) = JV<pi(x)V<f>f(x) < b\<fif >< <pf\S\4>i >< 4> t \a > 


( 2 . 8 ) 


Upon explicitly extracting the time dependence from the annihilation operators in 
m and (fkRl) by writing a k —> a k e~ luJkTi and 6 k —> 6 k e luJkT f , and using the functional 
integral representation of < 0/|5|0j >, one can write 


S(b*, a) = J V(f){x) exp iS[4>] + Rj(0i, a) + 5/(0/, b* 


(2.9) 


where the integral over the boundary configurations 0j and 0/ and the integral over 
all paths interpolating between these configurations have been combined into a single 
path integral, and the boundary terms are given by 
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Bi(<f>i,a) = / drk < --a k a_ k e 2 “ k * - - u; k 0j(k) </>;(-k) + '</>;(k) 


BMr.n = Jd 3 k{-\ kP- k ^ T ’ - ^ k 0,(k)^(-k) + Va^iie^^t-k) 


To proceed, the sums over initial and final states in (2T) are replaced by 


^2 ~ > / 'Da^Vcik exp 


d A k a£a k 


E 

/ 

and unity, in the form 


Vb^Vbk exp 


- / P£;6*6 k 


(2.11a) 

(2.11b) 


Vc^Vc^ exp 


— I d 3 /cc k c k 


Ic >< c| = 1 , 


( 2 . 12 ) 


is inserted between S' and PePn, giving 


a(E, N ) = J V[a , 6, c, e] exp [—6*6 — a*a — c*c — e*e] (2.13) 

5 ( 6 *, c)S(6*, e)* < c | PePjv | a > < a | P^P/v | e > . 

We are using an obvious short-hand notation for the integration measure, and in¬ 
tegrals over momenta are implied. The kernels of the projection operators take the 
form 


< 6|Pg|a > = J exp 

< b\P N \a > — dr) exp 


-iE£ + J d A ke iu ^ 6*a k 
—iNrj + J d A ke lv b^a^ 


(2.14a) 

(2.14b) 


from which it follows that 


< 6|PeP/v|o > = J dJ^drj exp *- iE£ — iNr] + J d A k e luJ ^ +vn 6 k a k 


(2.15) 


After substituting (|2.9|) and (|2.15|) into ( [2.131 ), and then changing variables via 
a —> exp[—zu;£ — i7]\ a and a* —> exp[— ioj^' — irf] a*, the c-integral may be performed 
to obtain 


(2.10a) 

(2.10b) 
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= exp [Bi(<f>i,a)\ , 


(2.16) 


J Vc*Vc exp 


c*c + c*a + c ) 


with a similar expression for the e-integration. Finally, after collecting terms and 
redefining £ + £' —> £ and 7] + r/ —> rj, the cross section becomes 

a(E,N)= j V<t){x)V(t)\x)ValVa v DblVb^dr 1 di exp [IT] , (2.17a) 

where 

IT = -*££ - iNr] - J d 3 k {fe*6 k + a£a k e" iAk } (2.17b) 

+ iS[<!>] - iS[<j>'} + Biifa, a) + B f {4> f , 6*) + 5,(0', a)* + Bf^, b*)* , 


with A k = a; k £ + rj. The functional 5[0(a;)] is the action, and the boundary terms at 
the initial and final times T) and Tj are given by ( p.lOj ). 

To display the semiclassical nature of the cross section, it is convenient to express 
the exponential factor W in terms of the rescaled held (j) = g(f>, the rescaled mode 
amplitudes a k = ga k and 6 k = gb k , and the rescaled energy e and particle num¬ 
ber v defined in ( |2.2| ). The action S'[</>], which is related to the unsealed action by 
S[0] = S[(ft\/g 2 , is ^-independent and we can thus write W(E,N ) = F(e, v)/g 2 , with 
the function F being independent of the coupling constant. Hence, for small values of 
g, it is a good approximation to simply saturate the integrals by classical saddle-point 
solutions, from which we obtain 


a(E, N ) = exp 



F(e,v) , 


(2.18) 


where F is determined by evaluating Q2.17b|) on the classical solution. In what follows 
we shall work only with the rescaled quantities in which the g-dependence has been 
factored out, but for notational simplicity we will use the unsealed notation and drop 
the tilde over the associated quantity. 
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In looking for saddle-points of (|2.17b| ) we must distinguish between two cases. 
There may be solutions which correspond to classically allowed evolution, in which 
case the fields and the action will be real and the parameter A k zero. As shown 


below, this implies that the function F of (|2.18|) will be zero, signalling the absence of 


suppression, and if we can find such classical solutions with small u, this furthermore 


indicates that the two-particle rates are likewise unsuppressed. Classically allowed 


evolution which changes the topology of the fields must perforce occur at an energy 
above the sphaleron barrier, but E > E sph is per se not a sufficient condition for 
the existence of classically allowed solutions with a given particle number in the 


initial state, and in a later section we shall return to the problem of finding topology 


changing solutions with low incident particle number. 

Alternatively, there may be solutions which correspond to classically forbid¬ 
den processes, in which case the saddle-points for 0 and 0', 5S[(p\/S(j) = 0 and 
8 S[(j)']/8<^' = 0, may in fact have imaginary components, while the saddle-point values 
of a and a* need not be complex conjugates. Obtaining these complexified saddle- 
points is much more involved than finding classically allowed solutions passing over 
the sphaleron barrier, and so we devote the remainder of this section to explicating 
some of the details of the procedure, with a special emphasis on boundary conditions. 

Extremizing (|2.17b|) with respect to modes Gk and a£ yields 


Gk e"* Ak + G_ k e~ 2iujTi - \[%j 00k) e~ iujTi = 0 


(2.19a) 


a k e"* Ak + a*_ k e 2iujT ‘ - >/2a;0j(-k) e iujTi = 0 , 


(2.19b) 


which may be solved to give the saddle-points 



(2.20b) 


(2.20a) 
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As previously noted, in general these solutions are not complex conjugates, and hence 
we use the bar notation for the latter. The expression Q2.20|) relates the a-mode 
amplitudes to the initial saddle-point values of the incident fields, which in turn are 
constrained by 


— k) — u4>i( k) + V2uj a_k e lwTi = 0 
i<j>i( k) - u;0'(k) + V2u al.e lujTi = 0 , 


(2.21a) 

(2.21b) 


and with the use of (|2.20|) , one can write this expression as a boundary condition 
involving only the incident fields, 


k) + u;00k) 



i0'(k) +o;0'(k) 
i0'(k) k) 


(2.22a) 

(2.22b) 


The parameter Ak itself is determined by saddle-point equations, and when it van¬ 
ishes note that 0'(k) = 0*(k), and that ak = (cnk/2) 1//2 0;(—k) e^ kTi and dk = 
(cn k /2) 1/2 00k) e~ luJkTi are complex conjugates. This case therefore corresponds to 
a classically allowed process above the sphaleron barrier. 

To obtain final-state boundary conditions, one extremes (|2.17b| ) with respect to 
the 6-modes and the final-state fields 0/ and 0)-, which gives 


6k + 6_ k e~ 2koTf - k) e~ iujTf = 0 (2.23a) 

6 k + 60 k e 2iuTf - VaJ 0/(-k) e iuTf = 0 , (2.23b) 

and 

k) —u(j)f(k) + V2ujhl i e luTs = 0 (2.24a) 

- k) - c^'/k) + 6_ k e~ iwTf = 0 . (2.24b) 


Together, (|2.23|) and (|2.24|) imply that the final-state fields and their respective time 
derivatives agree, 0j(k) = 0/(k) and 0y(k) = 0/(k). Thus, as the saddle-points 0 
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and (j)' satisfy the same classical equations at intermediate times, they also agree at 
these times, and hence we are really dealing with a single solution <p(x). At first sight 
this seems inconsistent with (|2.22| ) for nonzero A k . However, the general complex 
saddle-point solution is nonanalytic, and 0, and 0' in (|2.22|) are to be thought of as 
lying on separate sheets in the complex-/: plane (to emphasize this, we will not remove 
the prime from 0'). 

The value of A k is determined by the saddle-points of £ and /?, which from (|2.17b|) 
are related to the energy and particle number by 

e = y a£a k e _lAk (2.25a) 

v = J d 3 k <a k e~ iAk . (2.25b) 


The saddle-point of £ may be made pure imaginary by 
and the real part of the //-saddle-point is typically small 


a suitable time-translation, 
TT| so we can write 


i = iT 


(2.26a) 


r] — id . 


(2.26b) 


The parameter T can be removed from the boundary conditions by choosing the 
complex time contours of Fig. 1. Since the fields become linear in the distant past we 
can write 



FIG. 1. Complex-time contours 
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<f >{ k ) = 

0'(k) = 


■ 

1 

■ 


he 




+ 9k e 


ILJkT 




he 


+ 9ie u 


on line AB 
on line A'B' 


(2.27a) 

(2.27b) 


as r = Ret —> —oo. The boundary conditions ( |2.22| ) will be applied along line AB , 
keeping in mind that fa and fa r lie on different sheets in this region. In this asymptotic 
linear domain, however, the fields are analytic on their respective sheets, and hence 


0'(k) = 


1 




h e u}T ~ iuJkT + g' k e ~ ujT + iuJ k-r 


on line AB 


(2.28) 


as r —» —oo. This relation, with ( |2.27a| ) and boundary conditions (|2.22| ), gives the 
restriction 


h = e 6 h (2.29a) 

»k = e“V, (2.29b) 

and the degrees of freedom associated with the field along A'B' have been eliminated. 


A number of simplifications now occur. The energy and particle number may be 
written 

e = f d 3 k l n k h 5 -_ k (2.30a) 

v = Jd 3 khg- k . (2.30b) 

Upon taking the limits T % —> —oo and Tf —> oo the boundary terms become 

Bi(fa,a) = B*(fa,a) =\j d 3 kf^g^ (2.31a) 

B f (fa,b*) = B}(fa f ,b*) = l -j d 3 kb* k b k , (2.31b) 

and thus the exponential factor on the solution takes the form 

F = eT + v6 + iS[fa - iS[fa} , (2.32) 
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where the (rescaled g-independent) actions 5'[0] and S'[</>'] are evaluated along the 
upper contour ABCD on the first and second sheets respectively. 

By virtue of a symmetry akin to CPT , these expressions simplify considerably if 
the saddle-point is unique. Since the coefficients of the held equations 8S[(f>]/6(f> = 0 
are real, given a solution </>(x, t), one can form a new solution <l>(x, t) = 0(x, t*)*. 
Uniqueness then implies the conjugation symmetry 0(x, t) = 0(x, t*)*, and hence 
h = g *_ k and g k = /* k . This can be used to express (|2.29|) as 

9 k = e e fl k , (2.33) 


from which it follows that the energy and particle number take the form 


t = e - j (BkuJ'k ./k/ k 

u = e~ e ( d 3 kf*h . 


(2.34a) 

(2.34b) 


Expression 


may also be used to rewrite the boundary conditions (|2.27|) in 


the rather convenient form 


0(k) = 

= 


v2oJk ^ 
1 


he 




e e* iuj k r 


+ e ti he 


e° he^ + he 


* c “k7- 


on line AB 
on line A'B' 


(2.35a) 

(2.35b) 


\/2<Uk 

as r —> —cx). Note that the conjugation symmetry implies that the solution is real 
along the entire real-Minkowski axis, and we shall thus impose the additional bound¬ 
ary condition lm0(x, t = 0) = 0. In general, however, the solution becomes complex 
along the time-contours ABC and A'B'C , and the consistency of (|2.35a| ) and ( |2.35b|) 
requires that (j)(x ) also possesses singularities between these contours. For the case 
in which 6 = 0 , the held becomes real along AB and A'B' in the infinite past, and 
hence it remains real along the entire upper and lower contours. When the solution 
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is also real along the imaginary-time axis, as is the case for periodic instantons [9], 
turning-point boundary conditions 0 = 0 are also satisfied at B, B' and C. 

Recall that in ( [2.32D , both S = 5'[0] and S' = 5 , [0 / ] are evaluated along the contour 
ABCD , albeit on different sheets in the complex-f plane. If the singularities of 0 only 
lie between AB and A'B', the action S' along ABCD is equal to the action along 
A'B'CD (staying on the same sheet, and assuming no contribution from the contour 
at infinity), ffence, the conjugation symmetry implies S' = S*, and thus 

F = eT+ u9-21mS(T,9) , (2.36) 


where the implicit T and 9 dependence in the action has been made explicit. Note 
that S(T, 9) can be obtained by integrating only along ABC , as the contribution from 
the Minkowski section CD is real and does not contribute to F. Furthermore, since 
T and 9 (or equivalently £ and rf) are determined by the saddle-point of F, we have 


e = 21m 


v = 21m 


8S_ 

or 

dS 


86 ’ 


(2.37a) 

(2.37b) 


an alternate expression for the energy and particle number that can be used as a 
consistency check. 

In this section we have presented the formalism of Rubakov, Son, and Tinyakov 
[Q in some detail. The case of a single real scalar held has been used for purposes 
of illustration only, and the method is easily extended to more complicated theories. 
In most instances, and in particular for the standard model, the relevant saddle- 
point solutions must be obtained computationally. The function F(e, u) can then 
be determined, and the cross section ( j2.18| ) calculated. Finding these complexified 
saddle-points, however, is a formidable numerical challenge, and in the next subsection 
we outline some of our progress towards this goal and we present a few initial results. 
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C. Some Initial Computational Results 


Together with Peter Tinyakov, we are presently engaged in the formidable nu¬ 
merical task of finding the Euclidean-Minkowski hybrid solutions and extracting the 
baryon number violating cross sections. Our approach in the classically forbidden 
regime is to map the constant F-contours by exploring the two parameter family of 
solutions determined by 9 and T. As this investigation has just begun, we shall only 
present some preliminary results and briefly discuss our future plans. The numerical 
approach in the classically allowed domain above the sphaleron barrier is quite dif¬ 
ferent and involves exploring the e-u plane using Monte Carlo sampling techniques, 
and we shall present the details of this separate investigation in the next section. 

As previously stated, we are considering the standard model with the Weinberg 
angle set to zero. The resulting spontaneously broken SU( 2) gauge theory possesses 
all the relevant physics while undergoing notable simplification. The action for the 
bosonic sector of this theory is 

S — j dx 4 j-^Tr + D/FTFT - A($ t $ - 1) 2 | , (2.38) 

where the indices run from 0 to 3 and where 


Fpv = d^A v - d v A tl — i[Afj,, A u ] 
D li $ = (d li -iA li )$ . 


(2.39) 

(2.40) 


We use the standard metric r/ = cliag(1, — 1, — 1, — 1), and have eliminated several 
constants by a suitable choice of units. We have also set g = 1, but when needed we 
shall restore the gauge coupling to its physical value of g = 0.652. For our numerical 
work we take A = 0.1, which corresponds to a Higgs mass of about Mh = 72 GeV. 

To yield a computationally manageable system, we work in the spherical Ansatz 
of Ref. |H| in which the gauge and Higgs holds are parameterized in terms of six real 
functions ao , eq , a , (3 , n , and v of r and t: 
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4>(x,t) = -a 0 (r,t)cr • x 

A(x, t) = - [ ai (r, t)cr • xi* + ^^(<7* - <r • xfj + 1 + ^ *) +++V] 
2 r r 

$(x, t) = [pt(r, t ) + iu(r, t)cr • x]f , 


(2.41a) 

(2.41b) 

(2.41c) 


where x is the unit three-vector in the radial direction and £ is an arbitrary two- 
component complex unit-vector. For the 4-dimensional helds to be regular at the 
origin, a 0 , a, a\ — a/r, (1 + (3)/r and v must vanish like some appropriate power of 
r as r —> 0. 

These spherical configurations reduce the system to an effective 1+1 dimensional 
theory on the spatial half-line. The action of the reduced system follows by inserting 
(|2.41|) into (|2.38D, and after some algebra one obtains [|T^] 


S 


47T 



+ D'YA.V + T 1 D'‘n‘D i ,n, 


~ (Ixl 2 - l) 2 - \(\x\ 2 + 1)1+ - Re(YV 2 ) 
-Ar 2 (+-l) 2 l , 


(2.42) 


where the indices now run from 0 to 1 and are raised and lowered with rj lw = 
diag(l, —1), and where 


ffiu d^a v d t ,O j^ 

(2.43a) 

x = a + i(5 

(2.43b) 

(j) = ii + iv 

(2.43c) 

D^x = (d^-ia^x 

(2.43d) 

i 

= (<9 M - - a^cf) . 

(2.43e) 


This is an effective 2-dimensional f/( 1) gauge theory spontaneously broken by two 
scalar helds. It possesses the same rich topological structure as the full 4-dimensional 
theory and provides an excellent testing ground for numerical exploration. 
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In the rest of this section, we examine spherically symmetry Euclidean-Minkowski 
solutions lying on the ABCD time-contour of Fig. 1, satisfying the aforementioned 
boundary conditions. In particular, the real gauge held becomes complex and 
obeys the boundary condition 

a^(k) = g^k e~ luJkT + e e g*^ k e lUkT on AB as r = Ret —> — oo , (2.44) 

where we have absorbed the factor involving a4. into the definition of the amplitudes. 
For each value of the space-time index g, we are thus searching for two independent 
real degrees of freedom (as the saddle-point solution a M is complex). The real and 
the imaginary components of the complex fields y and 0 may be treated in a similar 
manner, giving a total of four real degrees of freedom for each held. We are thus 
looking for two independent complex helds y an d X, with boundary conditions 

x(k)=f k e- i ^ T + e e h k e i ^ T (2.45) 

X(k) = h* k e~ iuJkT + e e ft e iuJkT on AB as r -> -00 , 

and for two complex helds 0 and 0 with similar boundary conditions. When 9 = 0, 
the solution becomes “real” along the entire ABCD contour, in the sense that a M is 
real, while y = y* and 0 = 0*. 

In a future publication we plan to numerically find these solutions and to explore 
the behavior of the suppression function F(e, is) throughout much of the e-is plane. 
In this paper, however, we restrict ourselves to the periodic instantons of Ref. [9j] for 
which 6 = 0. 

As path EF of Fig. 2 illustrates, zero energy instantons are Euclidean solutions 
that interpolate between consecutive vacua. In contrast, periodic instantons are Eu¬ 
clidean solutions that interpolate between configurations of nonzero energy lying in 
adjacent topological sectors, as in path E'F' of Fig. 2. They have nonzero energy, 
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finite periods, and possess turning points (located at E' and F' in the Figure). Given 
a periodic instanton, we may choose the parameter T of Fig. 1 to coincide with the 
corresponding period, and hence we may take C and B to be the turning points E' 
and F' (with a time separation of iT/ 2). 



FIG. 2. The instanton ( EF ), the sphaleron (S), and the periodic instanton {E'F'). The hori¬ 
zontal axis represents the infinite dimensional field space, and the vertical axis marks the potential 
energy of a corresponding field configuration. 

We now wish to numerically find the periodic instanton along the imaginary time 
axis. These solutions are real along the Euclidean axis, in the sense that is real, 
while x = X* and (j) = (f)*- Moreover, since the time derivatives vanish at B and C, 
the periodic instanton remains real when continued both to the real axis and to the 
contour AB (consistent with the vanishing of 9). 

While finding these solutions is less challenging than obtaining the general saddle- 
points with nonzero 6, we must still resort to computational methods. However, before 
describing our numerical approach, it is useful to examine two instances in which F 
can be found analytically. Both the usual zero energy instanton and the sphaleron 
can be viewed as limiting cases of periodic instantons, and they are illustrated in 
Fig. 2. The instanton represents a vacuum-to-vacuum tunneling event, and as such 
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lies at the origin of the e-u plane, thus giving F(0, 0) = —2 Irn 5' inst = —167T 2 . The 
sphaleron, on the other hand, has both nonzero energy and particle number, e sp h and 
zz sp h respectively i. It is analytic over the entire complex-f plane, with a contribution 
to Inis' only along BC in Fig. 1, and hence F(e sp h, u sph ) = e sp h T — 2ImS' sph = 0. 

We now have the value of F at two key points in the e-u plane. The former gives 
the usual low energy’t Hooft suppression of the baryon number violation rates, while 
the latter yields unsuppressed rates at the sphaleron energy, albeit for very large 
incident particle number (as in a thermal plasma in the early universe). All other 
periodic instantons lie along a line connecting these two points, and F monotonically 
increases from — 167T 2 to zero as we traverse this line from the origin to the sphaleron. 
We now concentrate on finding these solutions and their corresponding values of F. 

It is convenient to work in the = 0 gauge, and to shift the zero of time so that 
turning points C and B are located at t — —iT/A and t — iT/ 4, respectively. That is 
to say, 

X(r, iT/ 4) = 0 (2.46a) 

4>(r,iT/ 4)=0 (2.46b) 

di(r,zT/4) = 0 , (2.46c) 

with vanishing time derivatives also at t — —iT/A. Since the fields merely retrace 
their steps after the turning points, it is sufficient to fold the periodic instanton only 
over the half-period from t = — iT/A and t = iT/A between consecutive turning points, 


1Technically, since the sphaleron is a static nonlinear configuration, it does not have a 
particle number; however, a gently perturbed sphaleron will decay into a state with a well 
defined particle number Z4 p h- It is the decaying sphaleron, or rather the time reversed 
solution, that we are actually speaking of here. To give a feeling for the numbers involved, 
when g = 1 and A = 0.1, the sphaleron energy is e sp h = 2.5447 and the asymptotic particle 
number is Z 4 p h = 1.7478 (in physical units with g = 0.065, E sp h ~ lOTeV and A sp h ~ 50). 
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which will considerably reduce our computational effort. In fact, we can be even more 
economical by exploiting an additional symmetry akin to time invariance, and it will 
then suffice to find the periodic instanton only over the quarter period, from t — 0 to 
t = iT/A. 

For any solution %(r, t), (ft(r,t), and ai(r,t) to the classical equations of motion, 
we may construct another solution given by 

X'(r,t) = ~X*(r, -t) (2.47a) 

ft{r,t) = —(ft*(r,—t) (2.47b) 

a[(r, t) = —aftr, —t) . (2.47c) 

This can be traced to a combination of parity symmetry (in the full 4-dimensional 
theory), (ft —» —(ft invariance, and simple time reflection t —» —t. In terms of the 
reduced theory, this is none other than the 2-dimensional time reversed solution 
(which should not be confused with 4-dimensional time reversal). 


A corresponding time reversed solution may be obtained from any given solu- 


tion by appropriately changing the initial conditions in a manner 

dictated by (|2.47). 

However, if the initial conditions take the form, 


x(r, 0) = imaginary y(r, 0) = real 

(2.48a) 

(ft(r, 0) = imaginary <ft(r, 0) = real 

(2.48b) 

o 

o 

(2.48c) 

(with no restriction on the time derivative of ai), then the solution will be invariant 

under time reversal: x' — X; ft — an d oft = cq, i.e. 


X(r,-t) = ~X*( r , t) 

(2.49a) 

(ft(r,-t) = — ft(r , t) 

(2.49b) 

ai(r, -1) = —a x (r , t) , 

(2.49c) 
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and in this case, we only need to look for solutions over the quarter period from t — 0 
to t = iT/4, satisfying the initial and final conditions (|2.46|) and (|2.48|). 


We must also choose boundary conditions at the origin and spatial inhnity in such 
a way as to ensure regularity of the corresponding 4-dimensional fields. For the y 
field, we must take y(0, t) = —i, and we may choose a gauge in which y(r, t) = i as 
r —> cx). Thus, as the initial ^-configuration is pure imaginary, it will always have a 
zero for some nonzero value of r. As will be discussed more thoroughly in the next 
section, the sphaleron is a spherical configuration with vanishing and with pure 
imaginary \ and 0; furthermore, like the initial configuration, its y-field possesses a 


zero at some radius r away from the origin. Hence, the initial configurations (|2.48|) 
will not lie far from the sphaleron, as indicated by the close proximity of S' and S in 
Fig. 2. 

To find the periodic instanton numerically, we place the system in a box of spatial 
extent L r and time extent L t , and then discretize the action ( |2.42| ) using the standard 
techniques of lattice gauge theory. The space-time grid has lattice sites at (iAr,jAt), 
where i = 0 ■ ■ ■ N r and j = 0 ■ ■ ■ N t . The fields y and </> become discrete variables 
defined on these sites, while cq is defined on the space-links and time-sites (in a 
gauge where a 0 does not vanish, it is defined on the space-sites and time-links). We 
take N r = 64, N t = 40, with dr = 0.05, and impose (|2.46|) and ( p.48|) on the lower 
and upper time slices j = 0 and j = N t respectively. The parameter At, which is 
taken between 0.02 to 0.04, controls the period of the periodic instanton through 
T = 4 A tN t . 


Starting from an initial guesses along the Euclidean axis satisfying the appropriate 
boundary conditions, we search for a minimum of the action using the method of 
conjugate gradients. Like the naive gradient descent, the conjugate gradient algorithm 
chooses its descent direction based upon the gradient. A new guess is then selected 
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lying further down the slope, and the algorithm is repeated. With each new iteration, 
the configuration finds itself closer and closer to the local minimum, and eventually 
one can approximate the extremum to within the desired tolerance. The advantage 
of the conjugate gradient method over a simple gradient descent, is that the former 
achieves a much more rapid convergence rate by judiciously shifting the direction of 
descent slightly away from the gradient. As this method is rather standard, we shall 
not describe it in any more detail. 

A straightforward application of the algorithm, however, yields little success. Pe¬ 
riodic instanton may be found in this manner, but unless one starts extremely close 
to a solution, one typically relaxes to the static sphaleron. The two turning points 
of the initial guess, E' and F' in Fig. 2, slide onto the sphaleron-like configuration 
S', which in turn eventually relaxes to the sphaleron itself. To avoid this, we add an 
additional term to the action that tends to pin the turning points, thereby halting 
the collapse into the sphaleron. We do this by minimizing an effective action of the 
form 


S eB = S + w t (V Nt -v o) 2 , (2.50) 

where V^ t is the potential energy on the final time slice, w t is a weighting factor, 
and u 0 is called the turning-point energy parameter. While the second term in (|2.50| ) 
renders collapse to the sphaleron energetically unfavorable, in general the minima of 
S'efr are not the solutions we seek. However, if we choose the parameter v 0 such that 
V Nt = v 0 , the minima of S e g and S coincide. Thus, by adjusting v 0 accordingly, we 
can find periodic instantons that do not collapse to the sphaleron. 

The effective action fl2.50| ) still typically fails to yield nontrivial periodic instan¬ 
tons. While the second term in Q2.50| ) pushes the gradient search away from the 
static sphaleron, there is still another unstable direction. Rather than converging to 
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periodic instanton solutions, most initial guesses shrink to zero size. This is related to 
the fact that the standard model, with a nonzero vacuum expectation value, breaks 
conformal invariance and strictly speaking does not support instanton solutions (i.e., 
they shrink to zero size). We can remedy the situation by adding another term to the 
action that pins the zero of y(r, 0), thereby halting the collapse of the configuration. 
We now consider the effective action 


SeS — S + W t (V)v t — Vo) 2 + V4ero [(1 — a )Xi ,0 + C%+l,o] 2 j 


(2.51) 


where w zero is another weighting factor. Just as before, we seek to minimize S e g, but 
in addition to adjusting v 0 so the second term of ( |2.51|) vanishes, we also vary a and i 
to give a vanishing third term (and hence the zero of y (r, 0) occurs at r 0 = (i + a)Ar). 

A conjugate gradient minimization of the effective action ( |2.51|) , coupled with the 
two parameter search over an d ct, is a very effective method for obtaining periodic 
instantons. In this paper we only have space to present a typical solution, shown in 
Fig. 3, and a comprehensive treatment of periodic instantons in the e-u plane must 
wait for a future publication. 



FIG. 3. Periodic Instanton: a full period of the %-field. The modulus of the field is represented 
by the height of the surface, while shades of gray code the phase. 
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III. THE CLASSICALLY ALLOWED DOMAIN 


We now examine the complimentary classically allowed regime above the 
sphaleron barrier, in which the solutions are purely real and propagate in Minkowski 
space-time. Finding these solutions is less computationally demanding than solving 
the tunneling problem of the previous section, while still yielding considerable infor¬ 
mation about baryon number violation. Spatial limitations prevent us from giving 
a full blown treatment of our numerical investigation, and the reader is referred to 
Ref. [§] for complete details. But the basic idea is that if a topology changing classical 
solution with small incident particle number could be found, this would be a strong 
indication that baryon number violation would be observable in high energy two- 
particle collisions. Conversely, if there are no small-multiplicity topology changing 
solutions, then it is unlikely that the rates become exponentially unsuppressed. 

This can be made more precise in the following manner. Because of energy dissi¬ 
pation, the system will asymptotically approach vacuum values and will consequently 
linearize in the past and future. Field evolution then becomes a superposition of 
normal mode oscillators with amplitudes a n , which allows us to define the asymptotic 
particle number v = Y \^ n \ 2 - Furthermore, since the fields approach vacuum values 
in the infinite past and future, the winding numbers of the asymptotic held config¬ 
urations are also well defined, and fermion number violation is given by the change 
in topology of these vacua [14|]. Because of the sphaleron barrier, classical solutions 
that change topology must have energy e greater than that of the sphaleron. The 
problem we would like to solve, then, is whether the incident particle number v of 
these solutions can be made arbitrarily small. That is to say, we wish to map the 
region of topology changing classical solutions in the e-v plane. 

We could easily parameterize incoming configurations in terms of small pertur¬ 
bations about a given vacuum, but it would be extremely difficult to choose the 
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parameters to ensure a subsequent change in winding number. This is because topol¬ 
ogy changing classical solutions must pass over the sphaleron barrier at some point 
in their evolution, which is extremely difficult to arrange by an appropriate choice of 
initial conditions. So computationally, we pursue a different strategy. We will evolve 
a configuration near the top of the sphaleron barrier until it linearizes and the particle 
number can be extracted. The time reversed solution, then, has a known incident 
particle number and will typically pass over the sphaleron barrier thereby changing 
topology. Of course we have no obvious control over the asymptotic particle number 
of the initial sphaleron-like configuration; however, by using suitable stochastic sam¬ 
pling techniques, we can systematically lower the particle number while ensuring a 
change of topology. This will allow us to explore the e-v plane and map the region 
of topology change, the lower boundary of which should tell us a great deal about 
baryon number violation in high energy collisions. 

A. Topological Transitions 


Let us now put some flesh on the bones of the above discussion. As in the previous 
section, we still consider the standard model with zero Weinberg-angle, defined by 
action (|2.38|). As before, the coupling constant has been set to unity, but when needed 


it will be restored to its physical value of g — 0.652. We also take A = 0.1, which 
corresponds to a Higgs mass of about M H = 72 GeV. Again we restrict ourselves to 
the spherical Ansatz ( 2.41 ), and examine the effective 1+1 dimensional (7(1) theory 


(|2A2l) • 


Before investigating classical solutions of this effective theory, it is useful to first 
explore its topological structure, which is very similar to that of the full 4-dimensional 
theory. Vacuum states of the effective 2-dimensional theory are characterized by 
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|xl = |0| = 1 and ix* 0 = —1 (as well as D^x — D^cj) = 0). The vacua then take the 
form 


^/x vac 



Xvac ^ & 


(3.1) 


vac 



where the gauge function Q = Q(r,t) is required to vanish at r = 0 to ensure reg¬ 
ularity of the 4-dimensional fields. Furthermore, like the full 4-dimensional theory, 
these vacua still possess nontrivial topological structure. Compactification of 3-space 
requires that Q(r, t) —> 2nn as r —> oo, in which case the winding number of such 
vacua in the cio = 0 gauge is simply the integer n. Note that as r varies from zero 
to infinity, x winds n times around the unit circle while 0 only winds by half that 


amount. 


Since the winding number is a topological invariant, a continuous path connecting 
two inequivalent vacua must at some point leave the manifold of vacuum configura¬ 
tions. Along this path there will be a configuration of maximal energy, and of all 
such maximal energy configurations there exists a unique one of minimal energy [7j. 
This configuration is called the sphaleron and may conveniently be parameterized by 


a s P h ( r ) = 0 
Xsph(r) = i[2f(r) - 1] 
0 sp h (r) = ih(r) , 


(3.2) 


where the profile functions / and h vanish at r = 0, and tend to unity as r —> oo and 
are otherwise determined by minimizing the energy functional. The sphaleron energy 
is approximately M w /g 2 ~ 10 TeV, or e = 47 t (2.54) for A = 0.1 in the units we are 
using (and it depends very weakly on the Higgs mass). 
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FIG. 4. The x and 4> fields for a vacuum-to-vacuum topology changing transition in a gauge 
inconsistent with compactified 3-space. The scalar fields are traced in the complex plane as the 
spatial coordinate spans the entire axis. Figs, (a) and (c) represent two inequivalent topological 
vacua while (b) is the sphaleron barrier separating them. 


While the form of the sphaleron given by ( fi.2|) , in which vanishes and x and 4> 
are pure imaginary, is convenient for numerical work, it does have a slight peculiarity. 
Recall that compactihcation of 3-space requires the gauge function U to approach an 
even multiple of 27T as r —> oo. It is possible to relax this restriction, and it will often 
be convenient to choose a gauge with U —> (2n+l)7r as r —> oo, in which case x vac —> i 
and 0 vac —*■ ±i. This is precisely the large-r boundary condition of the sphaleron, 
which illustrates that (|3.2|) is inconsistent with spatial compactihcation. There is of 
course nothing wrong with this, and a topological transition of unit winding number 
change in this gauge is illustrated in Fig. 4. Rather than x winding once around the 
unit circle, it instead winds over the left hemisphere before the transition and over 
the right after the transition. The total phase change is still 27 t, as it must be since 
this is a gauge invariant quantity. 


Throughout most of this paper we shall use a gauge consistent with ( |3.2|) in 
which space cannot be compactified. From a computational perspective, this will 
allow perturbations about the sphaleron to be more easily parameterized. There will, 
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however, be times when it is more convenient to impose spatial compactification, but 
we will always alert the reader to such a change of gauge. 


B. Classical Evolution 


So far we have primarily been considering topology changing sequences of config¬ 
urations, not necessarily solutions of the equations of motion. We now turn to the 
classical evolution of the system. We will consider solutions that linearize in the dis¬ 
tant past and future, and hence those that asymptote to specific topological vacua. 
This allows us to define the incident particle number, and it makes clear what is 
meant by the topology change of a classical solution (namely, the change in winding 
number of the asymptotic vacua). 

The field equations are coupled nonlinear particle differential equations and must 
be solved computationally on the lattice. But before we present our numerical pro¬ 
cedure, we first formulate the problem in the continuum. The equations of motion 
resulting from the action (|2.42|) are 


d^r 2 f, u ) = z [D„x*X ~ X*D vX ] + 7^ [D v <j>*<j> - <$>'D h 


(3.3a) 


D 2 + \(\ X \ 2 -l) + l 


X = 


(3.3b) 


D^r 2 D^ + — (| x 1 2 + 1) + 2Ar 2 


- 1 


<t> = i x4>* 


(3.3c) 


The v — 0 equation in ( ]3.3a| ) is not dynamical but is simply the Gauss’s law constraint. 

To solve these equations, we must supplement them with boundary conditions. 
The conditions at r = 0 can be derived by requiring the 4-dimensional fields to be 
regular at the origin. The behavior as r —> 0 must be 
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ao = a 0i ir + a 0 , 3 r 3 + ... 

(3.4a) 

CLi = Gp o + CLi^r 2 + . . . 

(3.4b) 

a = air + a3r 3 + ... 

(3.4c) 

(5 = -1 + for 2 + ... 

(3.4d) 

/i = /iO + + • • • 

(3-4e) 

v = Vir + z/ 3 r 3 + ... , 

(3.4f) 


where the coefficients of the r-expansion are undetermined functions of time. The 
r-behavior of the various terms are determined by the requirement that it has the 
appropriate power of r = ( x 2 + y 2 + z 2 ) 1//2 to render the 4-dimensional fields analytic 
in terms of x, y and z. For example, ao must be odd in r since A 0 is proportional to 
ao<r ■ x = (ao/r)cr • x. In terms of x and (j) the boundary conditions at r = 0 become 


ao(0, t)=0 (3.5a) 

y(0 ,t) =-i (3.5b) 

Re <9 r </>(0, t) = 0 (3.5c) 

Im0(O,t)=O. (3.5d) 


There is another r = 0 boundary condition which arises from the requirement that 
cti — a/r be regular as r —> 0. In the notation of (|3.4|) , this condition can be written as 
a i i0 = «i, and once imposed on initial configurations it remains satisfied at subsequent 
times because of Gauss’s law. 

We turn now to the large -r boundary conditions. Finite energy configurations 
must approach pure vacuum at spatial infinity, and we may choose a gauge in which 


a M (r, t) -> 0 

(3.6a) 

X(r, t ) -> i 

(3.6b) 

0(r, t) -> i 

(3.6c) 

29 



as r —> oo. This choice of gauge does not admit spatial compactification, but nonethe¬ 
less it is numerically convenient since it is consistent with the simple parameterization 
of the sphaleron ( |3.2| ) . At times we will choose a gauge consistent with spatial com¬ 
pactification in which %(r, t) —> —i and (p(r,t ) —► 1 as r —> oo, but unless otherwise 
specified we will take the large-r boundary conditions to be (|3.6|). 

The field equations ( |3.3|) , together with boundary conditions ( |3.5|) and ( fl.6| ). may 
now be used to evolve initial prohles and to investigate their subsequent topology 
change. The evolution is performed by discretizing the system using the methods of 
lattice gauge theory, in which we subdivide the r-axis into N equal intervals of length 
Ar with finite extent L = NAr (in our numerical simulations we take N = 2239 and 
Ar = 0.04). The field theoretic system then becomes finite and may be solved using 
standard numerical techniques. 

The fields x(r, t) and (p(r, t) become discrete variables Xi(£) and <pi(t) associated 
with the lattice sites r* = iAr where i = 0 ■ ■ ■ N. The continuum boundary condi¬ 
tions render the variables at the spatial end-points nondynamical, taking the values 
Xo = —i, Xn = i and (pjy = i (the value of </> 0 will be discussed momentarily). The 
time component of the gauge held ao(r,t) is also associated with the lattice sites and 
is represented by the variables do,i(£) with i = 0- • • N. We will usually work in the 
temporal gauge in which ao,i = 0, and we will not concern ourselves with this degree 
of freedom. 

The spatial components of the gauge field ai(r,t) become discrete variables asso¬ 
ciated with the oriented links of the lattice, and we represent them by a\^{t) = a;(f) 
located at positions r i+ 1 / 2 — (i + 1/2)Ar with i — 0 • • • N — 1. The covariant spatial 
derivatives become covariant finite difference operators that are also associated with 
the links, e.g. 

exp [—mi Ar/2] (p i+1 - (pi 


D r 


Ar 


i = 0- --N- 1 , 


(3.7) 
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where a* is short-hand notation for a^j. 

It is now straightforward to discretize the action (|2.42| ) in a manner that still 
possesses an exact local gauge invariance. But first, we need to state the restriction 
on (f>o(t) corresponding to the boundary conditions (|3.5c|) and (|3.5d|) . Since cq is real, 
we can write these boundary conditions in a covariant fashion by requiring the real 
part of D r <p and the imaginary part of (j) to vanish at r = 0. Using the discretized 


operator (|3.7|), we can then solve this boundary condition for </> 0 to obtain 


= Re [exp (—i a 0 Ar/2)(p 1 ] , 


(3.8) 


where a o is the value of Gq^ at i — 0 and should not be confused with the time-like 
vector held. This now allows us to eliminate 0o from the list of dynamical variables. 
Finally, the discretized Lagrangian becomes 


L = Att 


N -1 , r 2 


i =0 


i+ 1/2 


d 0 a,i - 


®0,i+l ,i 

A r 


exp(—f dj Ar )xi+i - Xi 
A r 2 


Ar 


N—l 


i=l 


+ 4v \ \(d 0 -ia 0>i )xi\ 2 + r*\(d 0 - -r* +l/2 


exp (-i at Ar/2)0 i+ i - & 


Ar 2 


(Ixil 2 + - l) 2 - Ar, 2 (l^l 2 - l) 2 j>Ar 


(3.9) 


2 [hn(exp(-f a 0 Ar/2)0i)] 2 
— 47T r 1/2 — , 

and the system may now be evolved using standard numerical techniques of ordinary 
differential equations. The Lagrangian ( |3.9| ) is actually of a Hamiltonian type with 
no dissipative terms, so it is convenient to use the leapfrog algorithm to perform the 
numerical integration. 

We do not have space to outline this well known computational procedure, so 
instead we simply state some of its more attractive features. First, the algorithm is 
second order accurate (i.e. the error from time discretization is of order (At) 3 in the 
individual steps and of order (At) 2 in an evolution of fixed length L = N Ar). Second, 
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energy is exactly conserved in the linear regime, a desirable feature when pulling out 
the particle number. And finally, the algorithm possesses an exact discretized-time 
invariance, which is important since we are interested in obtaining the time reversed 
solutions starting from perturbations about the sphaleron. Of course these last two 
properties hold exactly only up to round-off errors, which can be made quite small 
by using double precision arithmetic. 

C. The Initial Configuration: Perturbation About the Sphaleron 

We are now ready to continue our investigation into the connection between the 
incident particle number of a classical solution and subsequent topology change. We 
could proceed by parameterizing linear incoming configurations of known particle 
number, but it would be extremely difficult to arrange the classical trajectory to 
traverse the sphaleron barrier. If we failed to see topology change for a given initial 
configuration, we could never be sure whether it was simply forbidden in principle by 
the choice of incident particle number, or simply because the initial trajectory was 
pointed towards the wrong direction in field space. 

To alleviate this difficulty, we have chosen to evolve initial configurations at or near 
the moment of topology change, and when the linear regime is reached the particle 
number will be extracted in the manner explained shortly. The physical process of 
interest is then the time reversed solution that starts in the linear regime with known 
particle number and subsequently proceeds over the sphaleron barrier. Of course we 
must explicitly check whether topology change in fact occurs, but we have found that 
it usually does. Fig. 5 illustrates the numerical evolution of the x field for a typical 
topology changing solution obtained in this manner. The modulus of x is represented 
by the height of the surface, while the phase is color coded (but unfortunately we can 
only reproduce the figure in gray scale). We have reverted to a gauge where \n = — i 
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and (f>N — 1, consistent with spatial compactification, and in which the incoming 
state has no winding and the outgoing state has unit winding number. The topology 
change is represented by the persistent strip of 2tc phase change near the origin after 
the transition. 



FIG. 5. Topology changing transition: behavior of the \ field obtained the time reversal proce¬ 
dure described in the text. The various shades of gray code the phase of the complex field. The 
field starts as an excitation about the trivial vacuum, passes over the sphaleron and then emerges 
as an excitation about the vacuum of unit winding. Note the persistent strip of 2tt phase change 
near r = 0 after the wave bounces off the origin. 


We turn now to parameterizing initial configurations. For classical solutions that 
dissipate in the past and future, topology change (and hence baryon number violation) 


is characterized by zeros of the Higgs held |14|. For such topology changing solutions 
in the spherical Ansatz, the x held, which parameterizes the transverse gauge degrees 
of freedom, must also vanish at some point in its evolution. However, unless the 
transition proceeds directly through the sphaleron, the zeros of (j) and x need not 
occur simultaneously, and for convenience we shall choose to parameterize the initial 
configuration at the time when x vanishes for some nonzero r. Furthermore, we can 
exhaust the remaining gauge freedom by taking the initial x to he pure imaginary. We 
thus parameterize the initial conditions as an expansion in terms of some appropriate 
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complete set with coefficients c n , consistent only with the boundary conditions and 
the requirement that x be pure imaginary with a zero at some r > 0. 


We choose to parameterize initial conditions in terms of perturbations about the 
sphaleron given by linear combinations of spherical Bessel functions consistent with 
the small-r behavior (|3.4|). We only need the first three functions 


jo(x) = 
ji(x) = 


Sill X 


X 

sin x cos x 


x* 


x 


x° 


x 


J 2 (x) = ( —7 -) simc- — cos x , 


ar 


(3.10a) 

(3.10b) 

(3.10c) 


since jo(x) ~ 1, ji(x) ~ x and j 2 {x) ~ x 2 at small x. We also require the perturba¬ 
tions to vanish at r — L consistent with the large-r boundary conditions ( |3.6|) . We 
then parameterize perturbations about the sphaleron in terms of j nm (r) = j n (a nm r ) 
with n = 0,1, 2, where a nm with m = 1, 2, • • • are the zeros of j n (x). We are thus led 
to parameterize the initial conditions as 


^sph 

X(r, 0) = Xsph(r) c im hm(r) 

m= 1 

-^sph -^sph 

0(r, 0) = 0 sp h(r) + 5] C2m. jom ( r ) + * 5Z C 3m 3lm{ r ) 

m= 1 m= 1 

-^sph ^sph 

X(r, 0)=£ c 4m ji m (r ) +iJ2 c 5m j 2m (r) 

m= 1 m= 1 

-^sph ^sph 

0) = ^ 0)772 jomix') i ^ v ^7772 j 1?72 (^*) 

772—1 772—1 

^sph 

ai(r,0) = 51 C 8m 32m(r) , 

772—1 


(3.11a) 
(3.11b) 
(3.11c) 
(3.lid) 
(3.lie) 


where y sp h and </> sp h are the sphaleron prohles, and where the sum is cut off at 
N sp h < N. To avoid exciting short wave length modes corresponding to lattice arti¬ 
facts, we shall take lV S ph ~ 1V/50 (in our numerical work, lV S ph = 50 for N = 2239). 
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We have used continuum notation, but (|3.11| ) is to be thought of as defining x 


and 0 on the lattice sites 7y and oq on the links r i+ i / 2 . The time derivative of aq is to 
be determined by Gauss’s law. 

D. Normal Modes and Particle Number 

We are now in a position to discuss the manner in which the asymptotic particle 
number is to be extracted. Recall that once the system has reached the linear regime 
it can be represented as a superposition of normal modes, and the particle number 
can be defined as the sum of the squares of the normal mode amplitudes. Since we 
have put the system on a lattice, we should properly calculate these amplitudes using 
the exact normal modes of the discrete system. However, since our lattice is very 
dense (N = 2239 with Ar = 0.04), it suffices to project onto the normal modes of the 
corresponding continuum system of finite extent L = N Ar, the advantage being that 
we can solve for the continuum normal modes analytically. We have checked that this 
procedure agrees extremely well with projecting onto normal modes of the discrete 
system (obtained numerically), so for clarity we present only the continuum modes. 

ft is convenient to work in terms of the gauge invariant variables of Ref. [|T6| . We 
write the fields x and 0 in polar form, 

x = -t[l + y] e w (3.12) 


h 

1 H— 
r 


e lv , 


(3.13) 


where the variables y and h are gauge invariant. We can also define the gauge invariant 
angle 


Z = 0-2r ) , 


(3.14) 
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and in 1+1 dimensions we can define a gauge invariant quantity 0 through 


T fixv 2 , 


(3.15) 


where e 0 i = +1 and /i, v run over 0 and 1. Rather than working with the six gauge- 
variant degrees of freedom y, 0 and we use the four gauge invariant variables h, 
y, ip and f. 

We wish to find the equations of motion for small linearized fluctuations about 
the vacuum. In gauge invariant coordinates the vacuum takes the form h vac = y vac = 
VVac — ^vac — 0, and we thus need only work to linear order in the variables. From 
Ref. | Til the normal mode equations are 


d u &* + 4A )h = 0 


+ \ + %) y = 0 


fr > M 


du -0 - 


(IV^- 

1 + \r 2 


+ “n ^ = 0 


jr 2 d^ + c liV d , '‘ip 

1 + \r 2 




+ -£ = 0 . 


(3.16a) 

(3.16b) 

(3.16c) 

(3.16d) 


Equation (|3.16a|) corresponds to a pure Higgs excitation characterized by mass Mh = 
2\/A, while (|3.16b|) - ([3.16d|) correspond to three gauge modes of mass Mw = l/\/2. 0 
Note that there are four types of normal modes. The first two are easily obtained 
by solving the independent equations (|3.16a| ) and ( p.!6b| ), while the last two can be 
found by solving the coupled equations (|3.16c|) and ( |3.16d|) involving -0 and £. A 
solution in the linear regime can then be expanded as a combination of these four 
modes and the amplitudes a,k n extracted, where k = 1, 2, 3,4 specifies the mode type. 
The Higgs and gauge particle numbers are defined by 


*Upon restoring the factors of g and the Higgs vacuum expectation value v, these masses 
take the standard form Mh = v2Xv and My/ = (1/2 )gv. 
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^higgs / v |^ln| 

n =1 

^mode 

^gauge = ^ ^ | ^2 n | "b | ^3 n \ | n | J ? 

n =1 

(3.17) 

(3.18) 

with total particle number given by 


V ^higgs ^gauge • 

(3.19) 


To avoid counting lattice artifacts we take the ultraviolet cutoff on the mode sums to 
be N mode ~ N /5 to 77/10. 



FIG. 6. Decay of a small perturbation about the sphaleron: behavior of the particle number in 
the four modes as function of time for lattice parameters N = 2239, A r = 0.04 and N mode = 200 
with A = 0.1. The physical particle numbers are obtained by multiplying the asymptotic values in 
the graph by 47T /g 2 ~ 30, which gives -/Vhi ggs ~ 8 and -A gauge ~ 45, for a total physical particle 
number of N p hy S ~ 53. 

Space does not permit a detailed exposition of this procedure, and one should 
consult Ref. 0 for full details. Here we must be content with Fig. 6, which displays 
the behavior of the particle number in the four normal modes as a function of time. 
As an initial state we chose a typical perturbation about the sphaleron as described 
in the previous section, and we see that its evolution quickly linearizes and settles 
down into a definite asymptotic particle number. 
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E. Stochastic Sampling of Initial Configurations 


Recall that our computational strategy consists in evolving a configuration near 
the top of the sphaleron barrier until it linearizes, at which point the particle number 
is extracted and the time reversed solution is then used to generate the topology 
changing process of interest. We can regard the energy e and the asymptotic particle 
number v as functions of the parameters c n that specify the initial configuration, and 
by varying these coefficients we would like to explore the e-v plane and attempt to map 
the region of topology change. In particular, for a given energy e, we would like to find 
the minimum allowed particle number /y min (e) consistent with a change of topology. 
If this number can be made arbitrarily small, this would be a strong indication that 
baryon number violation would be observable in a two-particle collision. 

By randomly exploring the initial configuration space, parameterized by the co¬ 
efficients c n , we would stand little chance of making headway. Instead, we shall 
employ stochastic sampling techniques, which are ideal for tackling this type of multi¬ 
dimensional minimization. Our procedure will be to generate initial configurations 
weighted by W = exp(— F) with F = (3 e + fiu, and by adjusting the parameters (3 
and ft we can explore selected regions in the e-v plane. In particular, by increasing /i 
we can drive the system to lower and lower values of v for a given e. In our numerical 
work we typically take f3 between 50 and 1000 while v ranges between 1000 to 20000. 

To generate the desired distributions we have used a Metropolis Monte-Carlo 
algorithm. Starting from a definite configuration parameterized by c n , we perform 
an upgrade to c n —> c' n = c n + A c n where Ac n is Gaussian distributed with a mean 
of about 0.0008. We evolve the updated configuration until it linearizes and then 
calculate A F = f3 Ae + // Av. If the topology of the physically relevant time reversed 
solution does not change, then we discard the updated configuration. Otherwise we 
accept it with conditional probability p = Min[l, exp(—AF)], which is equivalent to 
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always accepting configurations that decrease F while accepting those that increase 
F with conditional probability exp(—A F). 



FIG. 7. Monte Carlo results with lattice parameters of N = 2239, A?’ = 0.04 (giving L=89.56), 
-Amode = 200 and N s p h = 50, and with a Higgs self-coupling of A = 0.1. The solid line marks 
the sphaleron energy e sp h = 47 t( 2.5426), below which no topology changing process can lie. The 
triangle represents the configuration from which we seeded our Monte Carlo search. To obtain 
quantities in physical units, multiply the numbers along the axes by 4ir/g 2 ~ 30. The energy 
axis extends from about 10 TeV to 15 TeV, while the particle number axis ranges from about 30 
particles to 60. 

We are now in a position to present our numerical results. Fig. 7 represents 300 
CPU hours and involves 30000 solutions (of which only 3000 are shown) obtained 
on the CM-5, a 64 node parallel supercomputer. We have chosen the lattice pa¬ 
rameters N = 2239, A r = 0.04, with ultraviolet cutoffs determined by 7V sph = 50 and 
N mo de = 200. The Higgs self-coupling was taken to be A = 0.1, which corresponds to 
a Higgs mass of M H = 72 GeV. 

We have managed to produce a marked decrease of about 40% in the minimum 
particle number // min (e), which is approximated by the lower boundary in the Fig. 7. 
Nowhere, however, in the explored energy range does v drop below 47T, or in physi¬ 
cal units the incident particle number N > 30 for energy E < 15 TeV (the outgoing 
particle number tends to be about 50 to 100). This is a far cry from two incoming par- 
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tides which would be necessary to argue that baryon number becomes unsuppressed 
in high energy collisions. 

The complex nature of the solution space can be illustrated by the break in popu¬ 
lation density between e/47r ~ 3 and e/47r ~ 3.4. In our first extended search we did 
not check whether topology change actually occurred, trusting that the time reversed 
solutions would continue over the sphaleron barrier. However, we later found an entire 
region between e/47r ~ 3 and e/Ai: ~ 3.4 in which the solutions never left the original 
topological sector. We excluded these points and restarted our search procedure near 
e/47r ~ 3. A small discontinuity in the lower boundary with slightly lower particle 
number was produced, but we have still managed to approximate z/ min (e) remarkably 
well. 

We can extract more information from the system by investigating the asymptotic 
spectral distribution |a*, n | 2 as a function of mode number n. Before we started the 
search, our seed configuration (represented by the triangle in Fig. 7) linearized into 
a distribution that was heavily peaked about a small mode number n p k ~ 50 (with 
An ~ 50), corresponding to a frequency of u; p k ~ nnp^/L ~ 0.1. After the search 
the solutions underwent a dramatic mode redistribution. The amplitudes |cifc n | 2 of 
the linear regime peaked at higher mode number, n p k ~ 75 — 100, with a much 
broader distribution (An ~ 200). Clearly our search procedure is very efficient in 
redistributing the mode population density. 

While v remains large throughout the energy range we have explored, it is in¬ 
teresting to note that n min (e) maintains a slow but steady decrease with no sign of 
leveling off. To obtain an indication of the possible behavior of v m - m (e) at higher ener¬ 
gies, we performed fits to our data using functional forms which incorporate expected 
analytical properties of the boundary of the domain of topology changing solutions. 
The fits gave a particle number N = 2 at energies in the range of 100 TeV to 450 TeV. 
Of course we must explore higher energies before drawing define conclusions, but this 
is at least suggestive that particle number might at some point become small. 
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IV. CONCLUSIONS 


We have reviewed in some detail the semiclassical method proposed by Rubakov, 
Son, and Tinyakov (RST) for bounding the exponential behavior of the two-particle 
baryon number violating cross section in the standard model. There are two distinct 
regimes, one in which the solutions that saturate the functional integral are real, 
corresponding to classically allowed processes above the sphaleron barrier, and an¬ 
other regime in which the saddle-point solutions are complex, representing Minkowski 
evolution followed by Euclidean tunneling under a barrier. In both cases, a topology 
changing solution of small incident particle number would be a signal that two-particle 
baryon number violating rates would be observable in high energy scattering experi¬ 
ments. 

Finding the aforementioned saddle-point solutions is a formidable numerical task, 
but we have nonetheless made considerable progress, and in these lectures we have 
presented some of our initial computational results. In the Minkowski regime above 
the sphaleron barrier, we evolve nonlinear configurations at the moment of topology 
change until the system linearizes, at which point the asymptotic particle number 
can be extracted. The time reversed solutions, which have known incident particle 
numbers, will typically undergo topology change, and our computational strategy is to 
stochastically search the space of such topology changing solutions weighted for small 
incoming particle number. We have found that our numerical algorithm is extremely 
efficient in sifting configurations of smaller and smaller incident particle number. 
Starting with a generic perturbation of the sphaleron, which decayed into about 50 
particles, we have managed to lower the particle number by 40%, to approximately 30 
particles, while still maintaining topology change. Even though this number is rather 
large, we have only covered a narrow energy range, from 10 TeV to 15 TeV, but still 
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it is noteworthy that over this domain the data show a slow but steady decrease in 
incident particle number as the energy increases. In an effort to increase the rate at 
which we can collect numerical data, we have continued our search on a lattice with 
half the number of sites. This new lattice is sufficiently dense to ensure adequate 
linearization, but with much less CPU time, and we soon hope to have results well 
beyond the energy range we have explored to date. 

Computational methods for finding the saddle-points in the classically forbid¬ 
den regime are more involved than the above stochastic procedure. These solutions, 
which can become complex along the Euclidean-time axis, satisfy rather complicated 
boundary conditions, and we are still developing a procedure robust enough to End 
the general RST saddle-point. The periodic instanton of Ref. 0 is a special case of 
these tunneling saddle-points, and it has the advantage that it remains real along the 
entire Euclidean axis and satisfies rather simple turning-point boundary conditions. 
In these lectures we have been content with presenting a computational procedure 
to solve for the periodic instanton based on conjugate gradient minimization, and we 
have presented a typical numerical solution. 
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